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O . Abstract 
O 

\q ' New insight into the correspondence between Quantum Chaos and Random 

on' 

Matrix Theory is gained by developing a semiclassical theory for the autocor- 
*^ ' relation function of spectral determinants. We study in particular the unitary 

§ ' operators which are the quantum versions of area preserving maps. The rel- 

O ■ evant Random Matrix ensembles are the Circular ensembles. The resulting 

^ ■ semiclassical expressions depend on the symmetry of the system with respect 

to time reversal, and on a classical parameter fi = trU — 1 where U is the 
classical 1-step evolution operator. For system without time reversal symme- 
try, we are able to reproduce the exact Random Matrix predictions in the 
limit fi — ► 0. For systems with time reversal symmetry we can reproduce only 
some of the features of Random Matrix Theory. For both classes we obtain 
the leading corrections in ju. The semiclassical theory for integrable systems 
is also developed, resulting in expressions which reproduce the theory for the 
Poissonian ensemble to leading order in the semiclassical limit. 
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I. INTRODUCTION 



One of the most remarkable achievements of quantum chaology, was the observation 
that spectra of systems which are chaotic in the classical limit, obey universal statistics 
which follow the predictions of Random Matrix Theory (RMT) M. In most cases, the above 
correspondence was established by comparing a few statistical functions, mostly two point 
densities, or functions derived from them like the A3 or the S 2 statistics. The nearest 
neighbor spacing distribution P(s) is another example. The study of n-points (n > 2) 
spectral statistics which would test the quantum chaos - RMT correspondence more closely, 
are much more difficult to implement in practice. 

In a previous paper || we proposed a different measure, which is based on the study 
of the statistical properties of the spectral determinant, sometimes referred to as the sec- 
ular function, the spectral ( function or the characteristic polynomial: It is the function 
which vanishes if and only if its argument belongs to the energy spectrum. For matrices of 
dimension N the characteristic polynomial is 

N N 

p H (x) = det(Ix - H) = nO - An(#))) = JXffy. (1) 

n=l 1=0 

The statistical properties of this function, can be expressed in terms of the statistics of either 
the eigenvalues A n or the coefficients a/. Since the two sets of variables are functionally 
related, they are statistically equivalent. In practice, however, one cannot check the full 
spectral distribution, and therefore it is advantageous to study statistical measures which 
are based on other accessible quantities. The measure which was the subject of || was the 
autocorrelation function 

C(0 = AT ( / + °° w(x)p H (x + £/2)p H (x - £/2)dx) (2) 

where w(x) is a positive function with a finite support and unit integral. M is a normalization 
constant so that C(0) = 1. Analytical expressions for C(£) were derived in |J and [0 for 
the standard RMT ensembles. These results can also be written down in terms of (\ai\ 2 ) H . 
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There are other advantage to the study of the statistics of the spectral determinant, 
which are especially important for establishing the connection between quantum chaos and 
RMT. The Gutzwiller trace formula, which gives the semiclassical theory for the spectral 
density, is not a proper function. The derivation of spectral statistics based on this theory, 
should therefore be augmented by additional assumptions or truncation procedures || . In 
contrast, the semiclassical expression for the spectral determinant involves a finite number 
of periodic orbits, and therefore it converges on the real energy axis. The semiclassical 
spectral determinant preserves another important property, namely, it is explicitly real for 
real energies ( f7j, |§, 0). Thus, the semiclassical study of the statistical properties of the 
spectral determinant can be based on a relatively solid starting point. Last but not least, 
the semiclassical spectral determinant shares many of its properties with the Riemann Siegel 
expression for the Riemann ( function on the critical line. The autocorrelation function for 
the Riemann ( is the subject of a recent study of Cheung and Keating [|1(|. The rigorous 



results obtained for the Riemann ( case, provide some support to the physically reasonable, 
(yet mathematically uncontrolled) approximations made in the semiclassical theory to be 
discussed here. 

The behavior of C(£) can be intuitively clarified, by considering two extreme cases: An 
equally spaced (infinitely rigid) spectrum produces a correlation function which is a strictly 
periodic function of £, 

C(0=cos7r£ (3) 

Where £ is measured in units of the mean level spacing. On the other hand, a Poissonian 
spectrum with N spectral points, yields a positive correlation which decays to zero on a 
scale which is proportional to y/N . Thus, the lack of correlation between the energy levels, 
induces a slowly decaying correlation function. The canonical random ensembles all display 
level repulsion which induce strong correlations. It is expected, therefore, that the "incipient 
crystalline character" |TT[] of the spectrum of these canonical random matrix ensembles will 
manifest itself by oscillatory, yet slowly decaying correlation functions. The more rigid the 



spectrum, the more marked and persistent will be the oscillations of C(£). 

In the present work we would like to concentrate on the spectra of Unitary operators, and 
use the autocorrelation function of their spectral determinant to study the correspondence 
between quantum chaos and RMT. The spectral statistics of Unitary operators appear nat- 
urally in quantum chaology, when one studies, e.g., Floquet operators corresponding to time 
periodic hamiltonians, |12| , fi~3H or scattering matrices |14|. The quantization of classical 



area preserving mappings also involve quantum unitary operators. |15[ In cases where the 
classical phase space is compact, the corresponding quantum Hilbert space is finite, and 
then one deals with unitary matrices of finite dimension. This is the case we shall study 
here. We shall provide a new derivation of the results obtained in p|, and add some new 
results. 

In the next chapter we shall define the correlation function for the cases of interest here, 
and quote the RMT results which were obtained recently by the Essen group |J. We shall 
then turn to the semiclassical theory and assuming the system to be classically integrable 
or chaotic, we shall establish the conditions under which correspondence with Poissonian or 
RMT can be achieved. 

II. SOME RESULTS FROM THE THEORY OF RANDOM MATRICES 

The spectrum of a. NxN unitary matrix S consists of N unimodular eigenvalues e ldl , 1 < 
I < N . It is convenient to write the characteristic polynomial so that it is real on the unit 
circle 

Z s (lo) = e^" ) det(J - e~ iuj S) (4) 
Where e* e = det(— S). The characteristic polynomial can be written down as 

N 

Z s (u) = e^- e ^*i^ 1 (5) 

1=0 

The unitarity of S leads to the relations 
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(6) 



The autocorrelation function reads now 
1 r 2n 



T N < 


a, 


2 


>p 


T N < 

2^n=0 - 




2 >P 



(7) 



where < ■ >p stands for the average with respect to the spectral measure of the circular 
ensembles of interest: the orthogonal (COE, (3=1), the unitary (CUE (3 = 2) or the 
Poissonian (j3 = 0) ensembles. The < \a{\ 2 >@ are the coefficients in the Fourier expansion 
of the correlation functions. In the sequel we shall develop the semiclassical theory for these 
coefficients. 

In the semiclassical derivation, we shall make use of the following identities: 



N 



det(I — xS) = a n x r 

n=0 



( ^X k \ 

exp r£^J 



(8) 



where s n = trS n . We define the generating function 



G p (x, y) = (det(J - xS) det(J - y&)) p = (exp (- £ -(x k s k + y k sl)j j 



(9) 



so that 



Qn+r. 



-G (x,y) 



(10) 



x,y=0 



n\m\ d n xd m y 

The most important message is that the building blocks of the generating function are the 
traces s n = trS n . 

The ensemble averages < |a n | 2 >p were calculated in for all values of (3: 

n(N — n) 



< a n > 



1 + 



N+l 



< aj > 



1. 



(11) 



For the Poissonian ensemble one gets || 



< \0"n\ >/3=0 = 



(12) 



It is convenient to introduce the scaled correlation length r\ = £ and in the limit of large 
N the sums in ([?]) can be approximated as integrals, to give 
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3 /sinvrr/ 1 d 2 sinvrrA sinvn? 
^=1(77) = - + — 7T^ ; Cp =2 {rj) = 13 

2 \ 717] 7T (777^ 7T77 y 1X1] 

The correlation function for the Poisson ensemble reads 

C^ =0 (t7) = cos^ — 77 « 1 - -(— 7?) 2 ~ exp -i- . (14) 



N V 2 v iV (/ J 1 \ iV 

The Poissonian correlation function decays much slower than the correlation functions of 
the other ensembles. 



III. THE SEMICLASSICAL THEORY 

The quantum unitary operator S which we consider, is assumed to be the quantum 
analogue of an area preserving map M. acting on a finite phase space domain of area A. (For 
the present discussion we shall confine our attention to maps with the twist property. The 
semiclassical treatment can be extended to the general case). The phase space coordinates 
are denoted by 7 = (q,p). The map can be defined in terms of a generating function which 
is the action q') 

d$(q,q r ) , = d$(g, q') _ 

dq ' dq' 



The explicit mapping function 7' = Ai(^y) is obtained by solving the implicit relations ([T3J). 
The semiclassical expression for the matrix elements of S in the q representation is 



(-Y 


~<9 2 $(g, q')' 


\2lriJ 


dqdq' 



1 

2 



(16) 



In the semiclassical limit, A, the dimension of the Hilbert space where S acts, is the 
integer part of 

We have seen above, that one can express the coefficients ai in terms of the traces 
s n = trS n . The semiclassical approximation for s n involves the periodic manifolds of the 
classical map. For hyperbolic maps [[TT)|], 

Sn = txS n « V ^ r (17) 



P n is the set of all primitive periodic orbits of M, with periods n t which are divisors of n, so 
that n = n t r. Orbits which are related by a discrete symmetry are counted once, and their 
multiplicity is denoted by gt- M t is the monodromy matrix and v t is the Maslov index. The 
action along the periodic orbit is 

nt 

= Z>(?j,?j+l) ( With 9n«+l = 9l)- ( 18 ) 
3=1 

For integrable maps, we use the phase space variables (/, 0) where / is the classical 
invariant. The domain of the mapping is / G [I m i n , Imax], £ [0, 2ir]. In this representation, 
the action takes the form <3>(0, 0') = $(0' — 0). The explicit form of the map is 

/' = / ; A0 = 0'-0 = F(/) (19) 

Where F(I) is the solution of the generating relation / = d ^j^f - ■ The expression for trS n is 
now 

n 
m=l 

Where the summation is carried over the periodic manifolds of period n and winding number 
m. They occur at values of / which satisfy F(I n , rn ) = 2tt^. 

Before we can proceed any further, we must clarify an essential point. In contrast with 
the RMT, the semiclassical theory deals with a single system, which does not depend on any 
random parameter which could represent an ensemble. However, averaging is mandatory in 
order to get a meaningful theory since the quantities we calculate fluctuate appreciably. We 
generate the "semiclassical ensemble" of s n , (with 1 < n < A/2) by considering the inverse 
Planck constant as a parameter, and different realizations of the ensemble are distinguished 
by the value of the parameter %~ l . We restrict h~ x to the interval I^T 1 — Tiq 1 ] < A . The 
mean value Hq 1 is large enough to justify the use of the semiclassical approximation. That 
is, for typical orbits ^o 1 ^* ~~ ®t'\ » 1- The interval A is taken to be small on the scale of 
Hq 1 , but sufficiently large so that A|$ t — $ t /| > 2n. In this way, the phases (mod27r) of the 
semiclassical expressions can be considered random. The averaging over the "semiclassical 
ensemble" is effected by 
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2tt 

nhF'{I ntm ) 



i[n$(A0=27rf)/7i-(n+i)f] 



(20) 



1 rh- 1 +A/2 

<A> h =- dh- l A(h- 1 ). (21) 

A Jn- 1 -A/2 

With this definition of the ensemble average, we get for both classically integrable and 
chaotic maps, 

< s n > h = 0. (22) 
The variance for the classically chaotic case reads, 



< Kl 2 >n« E i aTaTt i\/Tr\\ (23) 



teP„ 



det(/ - M t 



For integrable maps we get 



< \ s |2 > « Mima. Imrn) = ^ ^ 

2nn 

Thus, the Ti~ l averaging provides the well know diagonal (random phase) approximation. 
Due to the exponential proliferation of periodic orbits in hyperbolic maps, action difference 
become smaller as n increase, which is the reason why the diagonal approximation is not 
valid uniformly. We shall make use of the diagonal approximation in the restricted range 
n < A/2 m Hq 1 where the diagonal approximation is justified. For integrable maps, the 
variance of s n is independent of n. The result 4 < |s n | 2 1 implies that the spectral 

two-point correlation function for integrable systems is Poisson 0. 

At this point we would like to make a crucial observation: 
(4) The semiclassical ensemble of {s n } is an ensemble of independent random Gaussian 
variables. 

We shall show that this is approximately valid as long as the values of n are kept bellow 
the dimension A which is of order Kq . Consider the correlator < (s n ) k (s* n ) 1 >%. If n and m 
are relatively prime, the actions $ 4 which contribute to s n and s m are sufficiently different. 
Thus, all the terms in the product (s n ) k (s^ n ) 1 are oscillatory and will yield a vanishing 
result upon averaging. If n and m have a common divisor, j, choose k = m/j,l = n/j, 
and all the amplitudes which involve repetitions of the primitive orbits of length j will 
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contribute non oscillatory terms to the correlator. However, for hyperbolic maps, the number 
of periodic orbits which involve repetitions is exponentially smaller than the total number 
of periodic orbits, and the statistical independence of the variables s n and s m is ensured. 
The corresponding approximation is more difficult to justify for integrable maps because the 
proliferation of periodic manifolds is only algebraic. 

If we check all other correlators using the approximation that repetitions can be ne- 
glected, we find that the s n are Gaussian random variables. It should be noted that RMT 
implies that in the limit of large dimensions, (which coincides with the semiclassical limit) 
the traces are random Gaussian variables 0. This property, which is shared by the sta- 
tistical and the semiclassical ensembles, has far reaching consequences, and it constitutes a 
strong link between RMT and quantum chaos. 

Having established the statistical properties of the semiclassical ensemble we can calcu- 
late the generating function (§) for this ensemble. 



G h (x, y) = (det(I - xS) det(J - y&) 



(25) 



ex P -£7(^ + ^4) 



k=l 



exp J2(xy) 

\k=l 



k < \s k \ >h 

k 2 



so that 



1 d n 

< Clnttm >h = &n,m — 7 ~ Z%\V 1 

n\ o n v 



(26) 



v=0 



with 



Z h (v) = exp 1^2 < \ s k\ 2 >7 



\k=l 



k 2 



(27) 



The coefficients in the Taylor expansion Z%(v) = 1 + J2kLi can be obtained by taking 
successive derivatives of the two expressions for Z h (v) 



(28) 



We would like to emphasize the following important points 



• One can use (p8|) to calculate < |a/| 2 >%= A\ in the range 1 < I < A/2 exclusively. We 
have to impose this restriction because the derivation which leads to (|28| ) is not exact. The 
approximations destroy the delicate balance between the < |sfc| 2 >n which is necessary to 
maintain the symmetry < \a>i\ 2 >n=< \cla-i\ 2 >n which follows from unitarity. Because of 
the same reason, A\ ^ for I > A, which contradicts the basic ingredient of the theory, 
namely, that we deal with a characteristic polynomial of order A. 

• Because of the reasons listed above, 

C(0 + e~^ 2A Z h (e^) (29) 

• The input necessary to compute < \ai\ 2 >f l consists of all the < \sk\ 2 >n with k < I < A/2. 
Thus, the semiclassical ensemble of s n can be confined to the range n < A/2 where the 
arguments which support the conjecture can be most trusted. 

The semiclassical expressions for < \sk\ 2 >n, depend on the underlying classical dynamics 
(^),(^). We shall treat the various cases separately. 

A. Classical chaotic dynamics 

It is instructive to consider the Fredholm determinant for the classical evolution (Frobe- 
nius Peron) operator 

U(>yrf)=6(i-M(i)), (30) 

where 7 is a phase space point in the domain of Ai. A straight forward integration gives 

Un = trU n = V , , -r. (31) 

t ^Jdet(J-Mni 1 J 

Comparing this expression with (^3|) ; we may write 

< |S "' >SK «§„ |det(j'-Ag)l =< *"* >C ' ( "' ' <32) 
where for any phase space function ^(7), < p >d(n) is t ne approximate phase-space average 
of p obtained by sampling p on the set of periodic points of period n with the appropriate 
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weights. As n — > oo, < p > c u n ) approaches the true phase space average of p, since u n — > 1 



in this limit [17]. Defining 



g = hm < mt >d <"> (33) 
n 

we get 

< |s„| 2 >a~ ffnttf, (34) 

Thus, fl27|) can be written as 



Z^u) = exp ^ < |s fc | 2 > h 



\k=l ,v / 

« exp \22gu k — j (35) 

= exp (— gtrlog(X — vU)) 

= {det{I-vU)y 9 

= ((Ruelle(v)) 9 

Where (R ue u e ( z ) is the Ruelle ( function for the classical mapping M.. It is defined in terms 
of the Fredholm determinant of the classical evolution (Frobenius Peron) operator by 

(RueUe(z) = (det(I - ZU))- 1 . (36) 

The relation (|35D is one of the central results of the present work, and we shall come back 
to it after obtaining explicit expressions for < \a m \ 2 



Starting with the second line in (|35D we can easily derive the recursion relations 



< \a m \ 2 >h= — J2< \a m -k\ 2 >n u k (37) 
m 



m 



k=l 

2 



which should be solved with the initial condition < |ao| >n= 1- 

Let consider systems which are strongly mixing and for which all transients die out on a 
short time scale. In such cases, we may replace the u k by their asymptotic value uu — 1 for 
all k. The recursion relation (|3"?D reads now 
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< \a m \ 2 >h— ~y2< \ a m-k\ 2 >n ■ (38) 
m 



k=l 



This is solved by 



. I i2 ^ 



A 



(39) 



m + g — 

v f- 1 / 

as can be checked by direct substitution. This solution can be used only for m < A/2. In 
the range m > A/2 it should be augmented by the identity 

< h| 2 > h =< |a A _/| 2 > h (40) 

which is a direct consequence of unitarity. These results can be compared with the results 
of RMT (0). 

Systems without time reversal symmetry (TRS) have g — 1, and < |a m | 2 >n= 1. Thus, 
the semiclassical result coincides with the prediction of the theory for for the CUE, (see 
(0)), < \a m \ 2 > p=2 = 1 . 

In chaotic systems with TRS, g = 2 and 

| 1 +1 for 1 < I < A/2 

< M 2 >n ~ { ■ (4i) 

I 1 + A-l for A > / > A/2 
This expression does not reproduce the RMT result for the COE case ( ||) 

<M 2 > co , = i + 'xV'V^ (42) 

However, for large A, where the semiclassical approximation is justified, the semiclassical 
result agrees with the exact expression in a domain of / values of size \/A in the vicinity 
of the end points of the / interval, I = and I = A. The deterioration of the quality of 
the agreement between the semiclassical and the RMT expressions when TRS is imposed is 
typical, and it is an enigma in the field of quantum chaos 0. 

So far, we discussed systems for which all transients die out on a fast time scale, which 
was imposed by setting Uk = trU k = 1 for all k. This is possible only when one eigenvalue 
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of U is 1 and all the rest vanish. In generic systems, the spectrum is not degenerate in this 
extreme way. Rather, beside the eigenvalue 1 which corresponds to the conserved phase- 
space measure, the spectrum is in the interval [0, 1], and it accumulates at 0. The rate of 
decay of transients is determined by the magnitude of the eigenvalues of U which are less 
then 1. To get the leading correction due to the non vanishing eigenvalues of U, one can 
expand the recursion relation (|37|) to first order in // = U\ — 1. One obtains in this way 
recursion relations for the correction to < \a m \ 2 >hbar- They are particularly simple for the 
cases with g = 1,2 and the corrected coefficients are 

< k| 2 >hbar = 1 + A* for g = 1 (43) 

= l + l + 2fd for g = 2 

The symmetry < \ai\ 2 >hbar=< \dA-i\ 2 >hbar should be implemented for I > A/2. Recently, 
the Essen group studied numerically the variances of the coefficients of the characteristic 
polynomial for the quantum kicked top |J. They checked systems with and without TRS, 
and their numerical results show systematic deviations from the RMT predictions which 



are consistent with the expressions (f44|). The numerical results for the case without TRS is 
particularly convincing. 



B. Classical Integrable dynamics 

For integrable maps we have fl2"l|) 

< \s n \ 2 > a ~ A (44) 
The resulting recursion relations for the coefficients < \a m \ 2 >n are 

A m ^ I |2 v. 

. I |2 ^ 11 \ a m-k\ /-Atx\ 

<KI>»=-E 1 («) 

k=l 

We were not able to find a close form for the solution of this equation. However, to 
leading order < \a m \ 2 — which coincides with the leading term of the result for the 
Poisson ensemble (TT2]). As can be checked by direct evaluation, the deviation between the 
semiclassical and the exact values occurs already in the expression for < |a m= 2| 2 >n- 
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IV. CONCLUSIONS 



In this work we derived the semiclassical theory for the autocorrelation functions of 
spectral determinants, by obtaining explicit expressions for their Fourier coefficients. The 
semiclassical theory provided two basic ingredients - it supported the conjecture (4) that the 
semiclassical ensemble of s n form a random Gaussian ensemble, and, it provided the variances 

< |s n | 2 >n for the various cases under study. The comparison with the corresponding RMT 
results is not uniformly successful, and we would like to correlate the degree of success to the 
accuracy by which the semiclassical theory provides the two ingredients mentioned above. 

For chaotic systems without TRS, the assumption (4) is well founded, and the semiclas- 
sical expression for < |s n | 2 >n coincides with < |s n | 2 >p=2- Hence the semiclassical theory 
for the autocorrelation coefficients matches exactly the RMT result. 

For chaotic systems with TRS, the assumption is well founded, but the semiclassical 
expression for < |s n | 2 >^ agrees with < |s n | 2 >/3=i only in the low n domain. Hence the 
semiclassical and the RMT results agree to leading order only. 

For integrable system the assumption (4) is not so well founded, but the semiclassical 
expression for < |s n | 2 >% agrees with < |s n | 2 >p=o, resulting again in agreement only to 
leading order between the semiclassical and the RMT results. 

The last point which should be mentioned is the natural occurrence of the classical Ruelle 
( in the theory. This phenomenon was recently observed by several groups which studied 
the correspondence between quantum chaos and RMT using other statistical measures [|J, 
||. We would like to emphasize, whoever, that it is not correct to use (Rueiie directly in 
the quantum theory, (see (|29|) ). Rather, the quantum character of the problem must be 
introduced by truncating the coefficients at / = A/2, and imposing unitarity via the relation 

< | a; | 2 >=< |aA-;| 2 >• Remembering that A plays the role of the Heisenberg time in the 
present theory, we can interpret the above statement as the analogue of the "Riemann Siegel 
look-alike" || symmetry which is a basic ingredient in the semiclassical theory of spectral 
determinants. 
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